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Hydrogen can be a renewable energy carrier and is suggested to store renewable 
energy and mitigate carbon dioxide emissions. Subsurface storage of hydrogen in 
salt caverns, deep saline formations, and depleted oil/gas reservoirs would help to 
overcome imbalances between supply and demand of renewable energy. Hydrogen, 
however, is one of the most important electron donors for many subsurface 
microbial processes, including methanogenesis, sulfate reduction, and 
acetogenesis. These processes cause hydrogen loss and changes of reservoir 
properties during geological hydrogen storage operations. Here, we report the 
results of a typical halophilic sulfate-reducing bacterium growing in a microfluidic 
pore network saturated with hydrogen gas at 35 bar and 37°C. Test duration is 9 days. 
We observed a significant loss of Hz from microbial consumption after 2 days 
following injection into a microfluidic device. The consumption rate decreased 
over time as the microbial activity declined in the pore network. The 
consumption rate is influenced profoundly by the surface area of Hə bubbles and 
microbial activity. Microbial growth in the silicon pore network was observed to 
change the surface wettability from a water-wet to a neutral-wet state. Due to the 
coupling effect of Hə consumption by microbes and wettability alteration, the 
number of disconnected H2 bubbles in the pore network increased sharply over 
time. These results may have significant implications for hydrogen recovery and gas 
injectivity. First, pore-scale experimental results reveal the impacts of subsurface 
microbial growth on Ho in storage, which are useful to estimate rapidly the risk of 
microbial growth during subsurface Hə storage. Second, microvisual experiments 
provide critical observations of bubble-liquid interfacial area and reaction rate that 
are essential to the modeling that is needed to make long-term predictions. Third, 
results help us to improve the selection criteria for future storage sites. 


KEYWORDS 


hydrogen subsurface storage, microbial effects, microfluidic pore network, hydrogen 
consumption, contact angle, wettability 


1 Introduction 


Underground hydrogen storage (UHS) has been proposed as a reliable and safe technology 
to store large quantities of hydrogen (H3), which is produced from a surplus of renewable 
electrical energy (power-to-gas). This stored H, can subsequently compensate the seasonal 
fluctuations in the supply and demand of renewable energy (Carden and Paterson, 1979; 
Muhammed et al., 2022). Current subsurface storage projects in Europe often concentrated on 
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storing H, in artificially constructed salt caverns, e.g., the HyUnder 
project (Landinger et al., 2014) because salt cavern storage is suggested 
to be suitable to short-to medium-term energy demand fluctuations 
(Heinemann et al., 2021a). 

For long-term and large-scale H, storages, depleted oil and gas 
fields and saline aquifers are potential storage formations (Heinemann 
et al., 2018; Hemme and Van Berk, 2018), due to the large storage 
capability and well-known geological structure from former 
exploration (Zivar et al., 2021; Jafari Raad et al., 2022; Muhammed 
et al., 2022). The elevation of H, concentration following injection, 
stimulate the growth 
microorganisms in the This 
implications for gas storage and withdrawal efficiency including 
microbial H, consumption, gas composition changes, clogging and 


however, may of hydrogenotrophic 


subsurface. can have adverse 


corrosion (Dopffel et al., 2021; Thaysen et al., 2021). Several recent 
reviews discussed the potential microbial risks in subsurface H3 
storage and highlighted the importance of monitoring and controls 
of microbial activity in subsurface environments (Heinemann et al., 
2021b; Dopffel et al., 2021; Zivar et al., 2021). Thaysen et al. (Thaysen 
et al., 2021) estimated the risk for microbial growth in 42 depleted oil 
and gas fields and showed that microbial growth and H consumption 
will likely occur in low-salinity and low-temperature reservoirs. 

Three microbial metabolisms are associated with the highest risks 
for implementing UHS: methanogenesis, sulfate reduction, and 
acetogenesis. The stored H, can be microbially converted to CH4, 
H2S, and CH3;COOH respectively. The activity of different microbial 
processes depends on the availability of microorganisms and electron 
acceptors such as sulfate or carbon dioxide in the reservoir (Dopffel 
et al., 2021). Especially, sulfate-reducing bacteria are commonly found 
in hydrocarbon reservoirs, which prefer temperatures around 38°C 
and near-neutral pH conditions (Schwartz and Postgate, 1985), but 
can also be active at high salinity of 24% (Ollivier et al., 1991) and 
temperatures above 100°C (Jorgensen et al., 1992). Sulfate-reducing 
bacteria can use sulfate as an electron acceptor to oxidize H, and 
generate H,S, causing permanent H, loss and gas contamination 
(Ebigbo et al., 2013; Dopffel et al., 2021; Thaysen et al., 2021; Zivar 
et al., 2021). Due to their activity and growth, the high accumulation of 
bacterial cells can lead to the formation of biofilms and cause pore- 
clogging (microbial-induced clogging) near injection wells, resulting 
in dramatic decreases in gas injectivity (Eddaoui et al., 2021). A better 
understanding of microbial processes pertaining to subsurface H3 
storage is essential for precisely estimating the microbial risks and 
successful implementation at large-scale UHS in the future. 

To date there is little experimental data on mechanisms of 
microbial H, consumption (Harris et al., 2007) and the effect of 
microbial growth on gas storage efficiency is not fully understood 
(Thaysen et al, 2021). In this work, we used a silicon-wafer 
micromodel with a pore pattern from natural sandstone for direct 
observations of the microbial-induced sulfate reduction at 35 bar and 
37°C, representing the storage conditions of a shallow aquifer or a gas- 
water transition zone in a depleted gas field (Lysyy et al., 2022). A 
halophilic sulfate-reducing bacteria was cultured with H, gas for 
9 days. Analysis of time-lapsed image series enabled quantitative 
studies of variation in pore-scale H, saturation, contact angle 
evolution and disconnection with microbial growth. Our results 
add new experimental data to explain the microbial processes 
during UHS, that are essential for modeling efforts which are 
needed to make long-term predictions and improve the selection 
criteria for future storage sites. 
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FIGURE 1 

(A) The pore network (26.95 x 22.50 mm) was split into 16 columns 

and 18 rows for measurement of local biofilm/H2 saturation and 
porosity. Units are shown as red boxes, and the blue rectangles represent 
the inlet, middle and outlet areas. (B) Binary image of the repeating 

unit of the porous network (1.65 x 1.14 mm) with black grains and white 
pore space. (C) Simplified schematic of the experimental setup to show 
the indented flow into the pore network. A more detailed description of 
the experimental setup can be found in our previous study (Benali et al., 
2022) 


2 Materials and methods 
2.1 Bacterial strains and growth conditions 


Desulfohalobium retbaense DSM 5692, a halophilic sulfate- 
reducing strain, was used as the model bacterium in this study. It 
can grow at pH ranging from 5.5-8.0, optimum temperature of 37°C 
and optimum salinity of 12% (Ollivier et al., 1991). Under anaerobic 
conditions, Desulfohalobium retbaense can utilize H, as an electron 
donor and sulfate as electron acceptor producing H3S, for growth 
(Spring et al., 2010; Tinker et al., 2022). The organism was grown in 
the DSM 499 high salinity growth medium with some modifications: 
Base medium was: 1 g/L NH,Cl, 0.3 g/L K,HPO,, 0.3 g/L KH,PO,, 
20 g/L MgCl,-6H,0, 100 g/L NaCl, 2.7 g/L CaCl,-2H,0, 4 g/L KCI, 
3 g/L Na2SO,, 1 mL/L trace element SL-10, 0.5 mL/L Na-resazurin 
solution (0.1% w/v), 0.3 g/L Na2S-9 H20. 1.5 mL of a dense bacterial 
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solution was inoculated in 15 mL of base media amended with 80 mM 
sodium lactate, 40 mM sodium acetate, 0.3% yeast extract and 0.3% 
peptone solution at 37°C for 7 days under anaerobic conditions (N3 
headspace). This 7-day old culture was then used and injected into the 
micromodel. 


2.2 Experimental setup and microfluidic pore 
network 


The experimental setup consists of a high-pressure micromodel with 
a Zeiss microscope (Axio Zoom. V16, Zeiss) illuminating with a S ring 
cold-light (CL 9000 LED) source. This microscope setup was able to 
capture the whole pore network (total 121 separate images) with a 
resolution of 4.38 um/pixel and acquisition time of 73s. Pore space 
temperature was kept constant at 37°C + 0.5°C by circulating warm 
water through internal copper tubes in the chip holder. Experimental 
pressure was controlled at 35 bar by a high precision plunger pump 
(Quizix Q5000-10 K), and a back pressure regulator (EB1ZF1 Equilibar 
Zero Flow) connected to a pressurized 300mL N, cylinder. The 
experimental setup has a high-precision injection system with a 
constant dead-volume (10 uL before the pore network). The silicon 
micromodel capable of withstanding pressure up to 150 bar, was 
designed based on a thin section from a natural sandstone and the 
production procedures are found in (Buchgraber et al., 2012). The pore 
network constitutes of 36 repetitions of a unique pattern (Figure 1). The 
pore network has a porosity of 0.61 and pore volume of 11.1 uL. For image 
segmentation, the pore network was split into 288 units and units have a 
porosity range between 0.59 and 0.66. Three areas of interest (a unique 
pattern) have been identified at inlet, middle and outlet sites (blue 
rectangles in Figure 1A). The bottom silicon and top borosilicate glass 
have strongly hydrophilic surfaces due to the fabrication procedure where 
silicon is oxidized to SiO}. Two open channels connect port 1 to 2 and 
port 3 to 4 across the injection and production sides of the micromodel. A 
more detailed description of the experimental setup and micromodel 
properties is found elsewhere (Benali et al., 2022). 


2.3 Experimental procedure 
The experimental procedure consists of the following steps: 


1. Preparation: The pore space was flushed with ethanol, deionized 
(DI) water, and H20, (ACS reagent, 30 wt% solution in water, 
Thermo Scientific™) to clean all surfaces and waste from previous 
experiments, followed by > 100 pore volumes (PVs) of DI water. 

2. Pressurization: The microfluidic system was pressurized to the 
operating pressure of 35 bar by flowing the growth medium at 
50 L/min from port 1 to port 4. The pressure was controlled by the 
back pressure regulator at the outlet. 

3. Bacterial inoculation: To remove residual water in the pore 
network, approximately 40 PVs of growth medium were injected 
at 20 uL/min into the pore network. Then, 10 PVs of precultured 
bacterial solution were injected at 6 uL/min for inoculation and 
bacterial attachment to the pore space surfaces, followed by an 18 h 
of shut-in period. The period of bacterial inoculation was 
approximately 19 h. 

4. H, drainage: H, was injected at 5 uL/min in the pore space for 4 h. 
The volumetric flowrate corresponds to a Darcy velocity of 
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(A) Pore space biofilm growth during bacterial inoculation, H2 
drainage and microbial growth with Hə processes in the pore network 
over time at 37°C and 35 bar. The biofilm saturation (Sp) was measured 
rom high-resolution images of biofilms in the pore network over 
time (Supplementary Figure S1), with measurement uncertainty less than 
%. Sp equals the number of pixels of biofilms in the unit dividing with the 
number of pixels of pore space. Biofilms in the pore network were mainly 
ormed during the bacterial inoculation stage and the saturation of 
biofilm in the pore network (Sẹ) was only 8-1074. During the bacterial 
inoculation process, some biofilms grew in the inlet line to port 1 (see 
Figure 1C). The flow of H2 gas brought biofilms in the inlet line into the 
pore network, resulting in an increase of Sp (up to 1.3-1074). In the 
meantime, the shear stress from gas flow also caused the detachment of 
biofilms in the pore network. When microbial growth in H2 gas, cells 
dispersed from biofilms caused a reduction of Sp. Note that the plot 
starts from bacterial inoculation and ends at bacterial growth with H2 for 
9 days. The horizontal axis uses a logarithmic scale. (B) Distribution of 
biofilms in the pore network after Hə drainage (Day 1). Biofilms were 
mainly formed in the inlet regions (near port 1 and port 2) and partly 
blocked the pore throats near injection. 


10.15 mm/min. After gas breakthrough, 100 PVs of H, were 
continually injected into the pore network at the same 
injection rate. 


. Bacterial experiment: After drainage (i.e., 100 PVs of H, injected), 


all the tubing connected to the pore network were flushed with 
growth medium to remove H, and the microfluidic system was 
closed and incubated at 37°C for 7-9 days by flowing DI water 
The 
spatiotemporal distribution of H, gas in the pore network was 


through bypasses to maintain the system pressure. 
visualized by the high-resolution microscope system. The bacterial 


experiment started immediately after H, drainage (step 4. above), 
with a duration of 7-9 days. 
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FIGURE 3 


Plots of the dynamics of H2 saturation in the whole pore network, inlet, middle and outlet areas during bacterial (black solid lines) and sterilized (red dash 
lines) experiments. Note that all plots start from the end of H2 drainage, indicating that there is no loss from dissolution here. H2 saturation (Sg) was quantified 
from the segmented 2D images with the maximum uncertainty of 12% (i.e., maximum error = 0.12 x Sg). Comparing the sterilized and bacterial experiment for 
the whole pore network (A), the effect of bacterial consumption can be observed with more rapid decrease and lower S, for the bacterial relative to the 
sterilized experiment. The gas saturation in the inlet (B) has a linear decrease during bacterial experiment, and a 34.5% reduction in Sg from microbial 
consumption (compared with 3.9% for sterilized conditions). The middle area (C) has the similar trend with the pore network that microbial cells consumed 
approximately 14% of H2 in the first 2 days and only 6.9% in the next 7 days. Due to the low gas saturation in the outlet (D) from gas drainage, the gas phase was 
fully consumed in 2 days. Without microbial cells, the H2 loss in the inlet and middle areas was 3.9% for 5 days and no Hz gas was trapped in the outlet area. 
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A baseline, sterile experiment (i.e., without microbial cells) was 
performed with the experimental protocol detailed above. Each 
experiment was repeated twice with comparable results, and we 
discuss the main results using a single microbial experiment in 
Section 3 below. The repeated experiment demonstrates the same 
governing mechanisms for wettability change, microbial-induced H3 
loss and disconnection, although with higher initial H, saturation and 
lower biofilm saturation. The reader is referred to the Supplementary 
Material for more detail. 


3 Results and discussion 
3.1 Microbial-induced clogging 
Microbial growth in porous media very often leads to biofilm 


formation and microbial-induced pore-clogging (Liu et al., 2019), and 
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this was observed during bacterial inoculation and H, drainage 
processes in this work. Biofilms formed predominantly in the inlet 
area in less than 0.1% of the pore space after 18 h of growth. This is an 
expected result because D. retbaense is not a strong biofilm forming 
bacterium. We observed that some biofilms, however, blocked the flow 
of Hz gas (Supplementary Figure S1) and reduced gas injectivity 
occurred. For instance, during H drainage, the shear stress from 
gas flow pushed some biofilms formed in the injection line into the 
inlet area, and some biofilms detached from the pore network. There 
was, however, no new biofilm formation during microbial growth with 
H3 gas and the biofilm saturation (Sp) reduced by 24% compared to Sp 
after drainage and culturing in H, for 9days (Figure 2). Our 
hypothesis is that the biofilm reduction arises from a transition of 
bacterial lifestyle from biofilm-mode to planktonic-mode in a H,-rich 
environment (Jefferson, 2004; Chua et al., 2014; Liu et al., 2019). This 
hypothesis is corroborated by earlier work (Liu et al., 2019) where 
biofilm detachment occurred in high nutrient environments, but 
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Contours plots of H2 saturation (Sg) in the pore network during bacterial growth at 37 + 0.5°C and 35 bar. The inlet area near port 1 initially (Day 0) had a 


higher gas saturation (0.61), while only a few H2 bubbles saturated near port 4 (Sg = 3.5-10~*). After microbial growth with H2 gas for 1 day, the gas saturation in 
the pore network was decreased from 0.201—0.146, where the average S, reduction rate in each unit was 0.052 per day with a standard deviation of 0.0866. 
After Day 2, the reduction rate in each unit decreased to 0.009 per day and the H2 bubbles in the outlet area were fully consumed by microbial cells. From 

Day 4, the consumption rate decreased to near zero revealing that the microbial activities were dramatically reduced after 4 days of culturing in H2. To 


illustrate the Sg changes more clearly, contour lines in plots were smoothed with a parameter of 9.84-10~*. 


cannot be generalized for growth on H, based solely on the reported 
results: whether the lifestyle change is due to either reduced EPS 
formation or an active detachment process, and if whether cell 
numbers increased during this process, need to be investigated in 
future experiments. Our observations are only valid for the bacterial 
strain studied. Different microbial strains need to be investigated in 
the future to verify if this effect of H-induced biofilm detachment is 
common for many bacteria. 


3.2 Microbial-induced Hə loss 
In this study, the main microbial process with UHS implications is 


hydrogenotrophic sulfate reduction (Ebigbo et al., 2013; Hemme and 
Van Berk, 2018; Gregory et al., 2019): 
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SO; + 4H, + H* > HS” + 4H,0 


The sulfate-reducing bacterium (D. retbaense) used H, (aq) as a 
source of electrons to reduce sulfate and generate H3S. Figure 3 plots 
the changes of gas saturation in bacterial (solid lines) and sterilized 
(dash lines) experiments at inlet, middle and outlet areas (cf. Figure 1). 
Biofilm formation in the inlet area reduced gas injectivity (cf. Section 
3.1) and initial H, saturation (Sgo) for the inlet area (Sgo = 0.46) 
compared with the sterilized experiment (Sgo = 0.60). Note that all 
plots start from the end of H, drainage. After 2 days of microbial 
growth with H, gas, the H, saturation (S,) decreased by 29.4% due to 
the bacterial consumption, compared with 3.7% for the sterilized case 
(reduction in the sterilized case is believed to be a combined effect of 
system leakage and measurement uncertainty). In the following 
4 days, the reduction of Sg in the bacterial experiment was low 
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(5-107, and similar to the sterilized case), indicating less microbial 
activity. At the end of experiment (9 days), 32.9% of the initially stored 
H, gas was consumed by the microbial cells. The inlet area shows a 
constant bacterial consumption rate (reduction of S, over time) of 
0.018 per day (constant rate for entire experiment), which is lower 
compared to the middle (0.03) and outlet (0.025) area for the first 
2 days. On average, the reduction in H, saturation was one order of 
magnitude larger for the bacterial case (0.024) relative to the sterilized 
case (0.0029). Similar results were obtained from the repeated 
experiment (Supplementary Figure S2A) that the consumption rate 
in the first 2 days is much faster than the following days. Note that the 
initial H, saturation for the inlet, middle and outlet was not the same, 
and the differences of consumption rate were related to the 
distribution of bacterial solution and H, gas in the pore network. 

The spatial distribution of S, per unit during bacterial 
consumption (Figure 4) showed that the gas phase in the inlet area 
initially was connected, compared with predominantly disconnected 
H3 bubbles in the outlet area due to the effects of snap off (Singh et al., 
2017; Lysyy et al., 2022). The variation of disconnected H, bubbles in 
the pore network will be detailed in Section 3.4. The microbial 
consumption rate in the pore network was strongly influenced by 
the interfacial area between gas and brine. Smaller Hz bubbles 
improved the H, availability for microbial consumption due to i) 
increased H, dissolution and ii) larger contact area between H, gas and 
the microbial cells. We postulate that H, availability is a driving 
mechanism for the observed consumption efficiency and rate in the 
pore space. Rapid bacterial consumption (within 7 h) was observed in 
the outlet area, occupied with disconnected H bubbles with large 
contact area. The consumption rate is also determined by activity of 
the microbial cells, exemplified by a drop in consumption rate 
(constant S,) observed in whole pore network and middle area 
beyond 2 days. We attribute the observed rate reduction to changes 
of the microbial growing environment. We assume that both sulfate 
(electron acceptor) and acetate (carbon source) were still available at 
high concentration in the solution. The loss in activity might therefore 
be a result of other unfavorable conditions. Bacterial growth in a 
confined system will lead to formation of bioproducts and sometimes 
toxins accumulated in the aqueous phase. For example, increased H2S/ 
HS" concentrations from sulfate reduction can be toxic to cells, even to 
sulfate-reducing cells. This, together with a noticeable pH increase 
during sulfate reduction, will be unfavorable for continued bacterial 
growth (Liu et al., 2019). We were unable to measure H3S in the 
effluent due to the small volumes associated with microfluidic 
experiments. Hence, the mechanisms for the reduced activity need 
to be investigated further. Cells in the inlet area were less affected, most 
likely due to the presence of biofilms that are known to protect 
embedded cells from environmental stresses like pH change and 
toxins in the environment (Flemming and Wingender, 2010). 

Due to the low solubility of H, in brine, H, loss from diffusion and 
dissolution will be small compared to microbial-induced loss (see the 
Supplementary Material) in agreement with previous numerical 
estimations (Carden and Paterson, 1979; Anikeev et al., 2021). 
Microbial consumption, however, is a major risk for H, storage at 
industrial scale, and specifically for the microbial sulfate reduction in 
our work. As sulfate compounds are common in subsurface reservoirs, 
deriving from the aqueous dissolution of mineral phases like anhydrite 
(CaSO,4) (Hemme and Van Berk, 2018), H, injection can trigger 
hydrogenotrophic sulfate reduction due to accelerated activity of 
sulfate-reducing microorganisms. H> can be microbially converted 
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to H2S (aq), resulting in H3 loss (up to 32.9% in this study) and 
reservoir souring. Our results show that H, consumption rates vary 
over time due to changes in active microbial cells and mass transfer at 
the interface. Therefore, microbial induced H, loss pertaining to 
underground storage is determined by the reservoir conditions (T, 
p, pH), microbial cells (activity and density) and mass transfer at the 
interface. Previous theoretical studies have shown that consumption 
can vary from 0%-17% of the injected H3, depending on the aquifer 
conditions and the involved microbial mechanisms (Thaysen et al., 
2021). We measured H, consumption rates at high initial nutrient 
content, high initial cell numbers and optimal physicochemical 
conditions for the studied strain, hence, investigating a worst-case 
scenario for microbial H, loss for the oligotrophic subsurface. On the 
other hand, in the subsurface some effects will favor microbial 
consumption, for example buffering minerals (carbonates) which 
will stabilize pH leading to a prolonged microbial activity. Also, 
under certain flow conditions in the subsurface fresh nutrients and 
cells might be introduced into the storage area, resulting in a faster 
consumption rate in the late stage of gas storage compared to our 
confined system. 


3.3 Microbial-induced wettability alteration 


The drainage of H, into the pore network containing two solutions 
with the same salinity resulted in non-identical gas distribution 
(Figure 4; Supplementary Figure S3) as a result of microbial 
induced wettability alteration. The wettability changes were 
investigated by measuring the contact angle in a three-phase (gas- 
aqueous-grains) system, which was quantified with image analysis 
using open-source ImageJ software. The gas-liquid interfacial area 
reduced to near minimal values when the H, gas was exposed to the 
bacterial solution for 1 day (see images in Figure 5, Supplementary 
Figure $2B). Results indicate that the bacterial cells altered the surface 
wettability from an initial water-wet (41° contact angle) to a neutral- 
wet (average 96° contact angle), compared with an unchanged contact 
angle (28°) for the sterilized experiment. We assume that bioproduct 
and bacterial adhesion to the solid surfaces are the dominant 
mechanisms for the angle shift to increased 
hydrophobicity. Microbial cells growing within the water film 


contact 


around the H, gas phase can easily attach to the SiO, surface of 
the micromodel. The cell surface of microbes has wetting properties 
more neutral wet than the SiO, surface; thus, the wettability is altered 
depending on the nature of the bacterial cell’s membrane structure. 
Similarly, the increased pH due to microbial growth can also result in 
increased hydrophobicity, that can be explained by the protonation 
state of the polar molecules and, thus the number of charged surface 
species at the interface of H/aqueous/SiO, (Santha et al., 2019). 
Wettability is one of the key parameters that determines the 
storage capacity and gas recovery (Bai and Tahmasebi, 2022; Pan 
et al., 2022) for underground gas storage. With the increase of contact 
angles induced by microbes, the capillary entry pressure decreases 
resulting in more favorable drainage displacement efficiency. The 
contact angle near 90° also lead to minimal interfacial surface area 
and favorable gas recovery due to reduction in capillary pressure (Lin 
et al., 2019; Lin et al., 2021). This in turn leads to decreased, snapped- 
off, and disconnected gas phase and an increase in H, relative 
permeability during imbibition. Hence, a large percentage of gas 
can potentially be recovered from porous reservoir rock when all 
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FIGURE 5 


grains ranges from 24.4°-114.1". 
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(A) The average contact angle in the bacterial experiment increased from 41°-96’ within 20 h, while there was little variation in the sterilized experiment 
(kept at a constant of 33°). After 4 days, microbial cells changed the wettability of the SiOz surface to a neutral-wet state with the average contact angle of 100° 
and standard deviation of 11°. (B) Images of the increase of contact angle from Day 0-Day 1. The contact angle (8) among Hz gas, bacterial solution and solid 


other factors are equal. Accordingly, microbial growth induced 
minimal surface areas in the hydrogen-brine-SiO, system may 
enhance the gas storage efficiency of UHS operations. 


3.4 Microbial-induced bubble disconnection 


During H, drainage, gas-phase continuity can be broken by snap- 
off, generating isolated gas bubbles and trapping gas in the pore space 
(Lysyy et al., 2022). Due to the low H, drainage rate (the capillary 
number is 2.04-10™ and Pelect number equals 4.76) in our study, the 
initial saturation of isolated gas bubbles in the pore network was less 
than 0.01 after H, injection and mainly located in the outlet area 
(Figure 3). As discussed above, the microbial effects of bacterial 
consumption and the increase of contact angle can reallocate the 
residual H, phase in the pore network (see Figure 6). Initially, H3 
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preferentially occupied the large, connected pore clusters with large 
pore throats. After bacterial consumption for 1 day, most of the 
residual H, remained in the same pores (red) but some movements 
into the neighboring pores were also observed (yellow). As the 
capillary entry pressure was reduced with the increase of contact 
angle, more and more H, bubbles were able to enter the nearby pores 
with smaller pore throats. The disconnection of gas phase (S,j) in the 
pore network was measured by dividing the number of pixels of 
disconnected bubbles with the number of pixels of pore space. The plot 
in Figure 6A indicates that the disconnection of the H, phase in the 
pore network increased up to 72.8% for 5h, while the sterilized 
experiment increased only 22.6% and stayed at a constant value 
until the end of experiment. 

Bacterial consumption can break the continuity of the gas phase 
and generate isolated gas bubbles, which further increases the contact 
area of H) and microbial cells and then accelerated the gas 
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(A) Saturation of disconnected Hz bubbles (Sgi) over time with 
bacterial (black, solid) and sterilized (red, dashed) aqueous solutions. 
(B) Dynamics of H2 bubbles over 24 h: initially the Hə gas was distributed 
in large, connected pore clusters with large pore throats (grey + 

red), and the aqueous phase saturated the remaining pores (blue + 
yellow). The H2 bubble distribution changed over 24 h (red + yellow), 
where red represents the Hə bubbles that remained in the same pore 
clusters and yellow represents Hə bubbles with new pore occupancy. 
The aqueous phase distribution after 24 h (blue + grey) increased relative 
to initial conditions. The black represents silicon grains in the 
micromodel, and remain stationary during the experiment. We assume 
that the redistribution of gas after 24 h is partly due to contact angle 
changes and increased hydrophobicity. 


consumption rate. Therefore, a decrease of Sg; was observed after 7 has 
a result of the fast consumption of isolated gas bubbles. The same 
trend was observed in the repeating experiment (see Supplementary 
Figure S2C). Our pore-scale results reveal that microbial consumption 
and wettability alteration are the two main mechanisms that 
determine disconnection and entrapment of the gas phase in the 
pore network. The increase of disconnection of H3 gas will reduce gas 
mobility dramatically in the receiving reservoir (Alhosani et al., 2021). 
Furthermore, the blockage of gas bubbles at the entrance will increase 
the flow resistance (Jian et al., 2022), causing a low gas recovery. Most 
of the trapped H, bubbles in the high permeability zones will be 
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immobile in the subsequent injection cycles, hence, are difficult to 
recover. 

The microbial effects on disconnection of H, gas can be 
summarized as the following: the increase of contact angles 
induced by microbes can reduce the capillary entry pressure and 
decrease the snapped-off and disconnected gas phase during 
imbibition; the stored H, phase, however, can be consumed by 
microbial cells forming new disconnected H, bubbles in the pore. 


3.5 Implications for H2 underground storage 
in porous systems 


The aim of UHS is to store large amounts of H, gas from surplus 
production of renewable electricity to mitigate seasonal fluctuations in 
energy supply and demand (Amid et al., 2016; Heinemann et al., 2021b; 
Dopffel et al, 2021). Because H, is an import electron donor for 
microorganisms, assessment of microbial effects for UHS is an inevitable 
step to estimate the potential risks for gas injectivity, H, loss and recovery 
for large-scale Hp geo-storage operations (Anikeev et al., 2021). In our study, 
we investigated one of the major H consuming processes in the subsurface, 
bacterial sulfate reduction under high salinity conditions. Overall, four 
microbial effects were observed in our experiments: microbial induced- 
clogging, H, loss from bacterial consumption, wettability alteration, and 
increased residual trapping of the H, phase. 

Microbial-induced wettability alteration appears to be a positive 
effect of microbial growth in a H, storage reservoir: increased 
hydrophobicity resulted in the presence of minimal gas-liquid 
surface areas, and high recovery efficiency, due to reduction in 
capillary pressure and high gas relative permeability. In contrast, 
the other observed microbial effects are expected to affect adversely 
the storage efficiency. 

Biofilm formation in the near-well regions may reduce H; injectivity 
and change the subsurface transport properties. The risk of H, loss from 
dissolution and diffusion in the brine is limited due to low solubility, but 
permanent loss from microbial consumption must be considered for 
economic reasons, particular for storage in aquifers due to the abundance 
of sulfate-reducing microorganisms and sulfate in the formation water 
(Zivar et al., 2021). The production of hazardous gases, such as H2S would 
lead to contamination of the injected H}. The increase of disconnected H, 
bubbles in the pore network will also speed up the microbial consumption 
process due to increased interfacial area and increase the flow resistance, 
resulting in a low gas recovery. 

Whether all these microbial effects will be at play during H3 
storage will depend on field specific factors like the indigenous 
microbial community, cell numbers, available nutrients and sulfate. 
In two major H, storage field trials (Underground Sun Project, 
Austria, and HyChico, Argentina) microbial activity was reported 
as an important factor for H, loss, but without indications of active 
sulfate-reducing microbes because both reservoirs were low in sulfate 
and the major microbial activity was methanogenesis (H, + CO, = 
CH4) (Bauer and Austria, 2017). However, a recent study showed 
significant microbial hydrogen consumption with sulphate reduction 
in original brine from a gas reservoir (Dohrmann and Kriiger, 2023). 
The overall risk will depend on the field specific microbial 
communities and field characteristics. That means that further field 
trials need to show if the mechanisms described in our work will be 
observed in the field. Our results help in understanding and 
interpreting field test data and to avoid significant H, losses. 
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4 Conclusion 


Our pore-scale analysis comprehensive 


characterization of microbial effects of a typical halophilic sulfate- 


provided a 


reducing bacterium during geological H, storage operations within a 
porous network. We showed that the microbial activity during H3 
storage may have important implications for H, recovery, gas 
injectivity and reservoir properties. Considering the small scale and 
small volumes of our micromodels and the optimal culturing 
conditions, the estimation of H, loss rates from microbial 
consumption cannot directly be linked to large-scale operations to 
quantify the extent of the effects in a reservoir. Core-scale experiments 
and simulations are needed to further estimate the net positive or 
negative effect of microbial growth. As UHS has received a lot of 
attention as a technically feasible technology to store renewable energy 
and thereby mitigate climate change, the microbial risks need to be 
explored in more detail to avoid storage failure or even storage souring 
by H3S. Future laboratory research should not only focus on single 
strains but also use original field communities from potential storage 
locations. Furthermore, the influence of Hz on microbial living styles 
and also on rock properties in the underground storage site must be 
investigated further. 
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